Fluctuating pancake vortices revealed by dissipation of Josephson vortex lattice 
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In strongly anisotropic layered superconductors in tilted magnetic fields the Josephson vortex 
lattice coexists with the lattice of pancake vortices. Due to the interaction between them, the 
dissipation of the Josephson- vortex lattice occurs to be very sensitive to the presence of the pancake 
vortices. If the c-axis magnetic field is smaller then the corresponding lower critical field, the 
pancake stacks are not formed but the individual pancakes may exist in the fluctuational regime 
either near surface in large-size samples or in the central region for small-size mesas. We calculate 
the contribution of such fluctuating pancake vortices to the c-axis conductivity of the Josephson 
vortex lattice and compare the theoretical results with measurements on small mesas fabricated out 
of Bi2Sr2CaCu2 08+5 crystals. A fingerprint of fluctuating pancakes is characteristic exponential 
dependence of the c-axis conductivity observed experimentally. Our results provide strong evidence 
of the existence of the fluctuating pancakes and their influence on the Josephson-vortex-lattice 
dissipation. 



I. INTRODUCTION 

Vortex physics in layered superconductors is extremely 
rich and interesting. This is mostly related to the fact 
that the magnetic field parallel to the layers penetrates 
in the form of the Josephson vortices (JVs) [T] while the 
perpendicular field creates the stacks of the pancake vor- 
tices (PVs) |2j . In the tilted magnetic field the attractive 
interaction between JVs and PVs 3J leads to the cre- 
ation of many unusual vortex states: mixed chain-lattice 
states 0] , JVs decorated by PVs [S] , etc (see as a review 
Ref. [6]). The interaction between the PV and JV lattices 
is revealed also in their dynamic properties. Dynamics 
of the JV lattice in Bi2Sr2CaCu208+5 (Bi2212) mesas 
have been explored by several experimental groups |3- 
I12j . Due to the weak intrinsic dissipation, the friction 
force acting on the moving JV lattice is small. That is 
why even the presence of a small amount of the PV stacks 
strongly affects the dynamic of the JVs. The physical rea- 
son is that the moving JVs induce displacements of the 
PVs which have large viscous friction due to the normal 
cores. This strongly enhances the JV friction leading to 
decrease of the measured c-axis resistivity of the sample. 
Recently, this effects has been studied in details in Refs. 
[T51 and [m To our knowledge, this is the only known 
situation in which adding vortices to a superconductor 
reduces its resistivity. Therefore, studying the dynamic 
properties of the JV lattice, we may obtain the informa- 
tion about the presence of the PVs. 

In the present paper we use this tool to study the 



penetration of the PVs in superconducting mesas of the 
Bi2212 intrinsic Josephson junctions. If the magnetic 
field is applied exactly along the layers, then the c-axis 
resistivity is governed by the friction force acting on the 
moving JV lattice. When the perpendicular component 
of magnetic field is switched on, the PV stack can 
enter into the sample when reaches the value of the 
perpendicular low critical field ffg^ P^ 

The c-axis field penetration is rather special in the case 
of layered superconductors with weak Josephson coupling 
between the layers. Indeed, the field starts to penetrate 
in the form of the PVs created near the surface in the 
fluctuational regime. The energy of a single PV near the 
surface has been calculated in Ref. [1^] (see also [T7]) and 
above certain field (somewhat lower than iJ^i) it has a 
minimum at distance of the order of London penetration 
depth A from the surface. In small-size samples with 
lateral dimensions comparable with the London penetra- 
tion depth this energy minimum appears in the center 
of the sample. The concentration of the fiuctuation PVs 
is determined by this minimum energy, which, in turn, 
depends on the perpendicular magnetic field. 

The presence of these fluctuational PVs increases the 
friction of the JVs and then it may be directly monitored 
by the measurement of the c-axis resistivity. We obtain a 
very characteristic exponential dependence of the c-axis 
conductivity vs , reflecting the varying density of the 
thermally activated PVs near the surface or in the center 
of small samples. This dependence is clearly observed 
experimentally. These results demonstrate that the JVs 
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serve as a perfect tool to detect the presence of the PVs 
in the fluctuational regime. 

The behavior of the JV lattice has special features in 
the mesas with lateral sizes of few Josephson lengths used 
in our experiments. In this case bulk interactions fa- 
voring the triangular lattice compete with edge interac- 
tions favoring the rectangular lattice. Due to this com- 
petition, the JV lattice undergoes a series of structural 
phase transitions between the triangular and rectangu- 
lar configurations. [T51 [T^ The regions of rectangular lat- 
tice dominate at large fields. As in the rectangular lat- 
tice the in-plane current are absent, its interaction with 
PVs is very weak. Therefore the fluctuational PVs have 
strong influence only on dissipation of the triangular lat- 
tice. Due to commensurability effects, the c-axis resistiv- 
ity has pronounced oscillations as function of the mag- 
netic field. 1 10) The period of these oscillations changes 
from half flux quantum per junction at small fields to 
one flux quantum per junction at large fields |1HI12| and 
the crossover field between these regimes is proportional 
to mesa width. jTHI 

The article is organized as follows. In section |Tl] we 
describe theory of the effect. In particular, in subsec- 
tion [HA] we derive the general formalism describing the 
response of the statistical ensemble of PVs in the static 
potential (creating by interaction with the surface) to the 
moving periodic potential (due to the interaction with 
moving JV lattice). We provide the expression for the 
resulting dissipation. In subsection |IIB 1| we detail the 
potential of the PV interaction with surface in presence 
of the perpendicular field . We also consider potential 
energy of a PV in finite-size samples with simple geome- 
tries (strip and cylinder). The energy of the interaction 
between PV located near the surface and the dense lat- 
tice of the JV's is calculated in subsection III CI We ana- 
lyze the general expression for the PV contribution to the 
JV lattice conductivity in section |II D[ In section |III| we 
present the experimental data on the c-axis conductivity 
of Bi2Sr2CaCu208+5 mesas with different lateral sizes in 
tilted magnetic field and compare the results with theo- 
retical predictions. 




FIG. 1. (color online) These schematic figures illustrate mesa 
geometry and the interaction between the dense JV lattice 
and fluctuational PVs near the surface, (a) Right part shows 
outline of the structure fabricated from Bi2212 crystal to 
study c-axis transport through highlighted mesa region. Left 
part shows definitions of axes and mesa sizes, (b) Three- 
dimensional sketch of the JV lattice inside the junctions (front 
face) and the fluctuating PVs near the surface (top face), (c) 
The graylevel plot of the cosine of phase difference between 
two neighboring layers which have periodic modulation due 
to the JV lattice and perturbations due to randomly located 
PVs in two layers shown by open and closed circles, (d) The 
energy profile of PVs near the surface in the c-axis magnetic 
field close to the lower-critical field. The PV density is en- 
hanced near the potential minimum. This highlighted region 
near the minimum of the potential gives dominating contri- 
bution to dissipation of the JV lattice. 



ity V = uj/k. This motion generates an electric field 
Ez = Hxv/c. Therefore, we are interested in the behavior 
of the PVs in the two-dimensional potential 



U{r, t) = Uo{x) + U{x) cos{ky - ojt). 



(1) 



II. THEORY 

A. Fluctuation pancake vortices in static and 
moving periodic potentials: General consideration 



We start from a general consideration of a model in 
which the PVs are considered as an ensemble of parti- 
cles with concentration n in the external time-dependent 
potential U{v,t). The particle current j [20] and particle- 
conservation equation are given by 



We consider the geometry with the parallel-to-layers 
magnetic field directed along x-axis iJy — and there- 
fore, the perpendicular current (along c-axis ) will drive 
the JV lattice along y-axis, see Fig. [ija). The PV lo- 
cated at distance x will interact with the surface (yz 
plane) and with a current screening the perpendicular 
field Hz as well as with a moving JV lattice. The first 
two contributions are described by the static potential 
Uo{x). The moving periodic potentials of the JV lat- 
tice is U{x) cos{ky — ujt), corresponding to the veloc- 



j = -D I Vn 
dn 



VU 



dt 



VC/ 

? 

T 



0, 



(2) 
(3) 



where D is the diffusion constant. 

In the absence of the moving potential the particle cur- 
rent is naturally equal to zero and the equilibrium density 
riQ^x) oc exv(—Un(x)/T).^IT\ Introducing the perturba- 
tion for the density v{r, t) and current due to the moving 



3 



potential 



n{r,t) ^noix)[l + u{r,t)], (4a) 
J=-D(viy+^^no, (4b) 

we have the following equation for the small perturbation 

^no - DV {noVu) - (noVU^ . (5) 

Using the complex presentation 

I7(r, t) = U{x) Re {exp [i{ky - ujt)]} , 
iy{r,t) = Re{vi^{x) exp [i{ky - ujt)]} , 

we obtain equation for the complex amplitude v^jix) 

T 



no V / 



no 



(6) 



where rj = T/D \s the viscosity coefficient. 

Assuming that the frequency is small, we solve this 
equation by iterations. The first iteration corresponds to 
quasiequilibrium 



The energy dissipation in steady state can be evaluated 
as 



W = - f dx{i\IU)=T] f dx(^- 



V I ^ Jj.\' + \jy\\ 



no 



Substituting currents (9a) and (9b I, we finally obtain 



(10) 



w = 



2 / C/„,i„\ 



— n.exp 1-— 



dxriijix) 

pi IJ2 



K \{K+Tk^f T^k^ 



(11) 



This formula determines energy losses caused by arbi- 
trary moving periodic potential near the minimum of the 
static potential. We will use it to evaluate the contribu- 
tion of the fluctuating PVs to the JV lattice conductivity. 
For this, in the following sections, we will obtain explicit 
expressions for the static and dynamic potentials. 



B. Energy of pancake vortex 

1. Energy profile of pancake vortex near surface in 
large-size samples 



The next iteration has to be found from equation 



D 



no V / 



U 

lUJ — . 

T 



(7) m 



Near the minimum of the potential the equilibrium den- 
sity can be represented as 



no(x) = nooexp 



Kx^ 



"-00 = exp (-[/niin/T) . 

We also assume that the typical scale of U{x) variation is 
larger than confining length of the static potential yjT jK 
so that we can use the expansion near the minimum, 
lJ(x) KiUo — Fx. In this case Eq. ([t]) has simple analyt- 
ical solution 



.(1) 



icuUo 



iLoFx 



Dk^T D{K + Tk^y 
which gives the oscillating particle-current amplitudes 



(8) 



K + Tk^ 
iiuU 
Tk 



"0, 



(9a) 
(9b) 



If the perpendicular magnetic field is slightly lower 
than H^^ , the PV stacks ( Abrikosov vortices) do not pen- 
etrate into the sample, but the PVs can exist near the 
surface at the fluctuation regime. The energy of the pan- 
cake vortex at the distance x from the surface is given by 



Uoix) 



.$2 



(47rA) 

5^0 (ggi-g 

47r 



2 



expi-- 



(47rA) 



47r 



expl-- 



(12) 



where s is the distance between superconducting layers 
and i/^^i w [$o/(47rA2)][ln(A/0 -hO.5] is the lower critical 
field. Above certain field this potential will have a min- 
imum at distance ~ A from the surface, see schematic 
Fig. [IJd). Taking derivative. 



dUo 
dx 



£$0 

47rA 



$0 A fx 



we find the condition for the minimum 

"""P^-") = 4^A^' 
with u = x/\. The minimum exists if 

e$o 



47rA2' 



(13) 
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Above this field the pancake energy at the minimum is 
determined by equations 



hr (Wmin) + 



1 



(14) 



exp (-■!/„ 



where Sp — s$o/ (47rA)^ is the PV energy scale. As Umin 
roughly behaves as ln(47rA^i/2/$o), the main dependence 
of C/min on is given by the first term, i.e., the minimum 
energy approximately linearly decreases with the field. 
For Bi2212 s « 1.56 nm and we estimate the slope of 
this dependence as dUmin/dH^ « —18.6 K/G. 



Expanding the potential (12) near the minimum, 
U{x) — Ui-nin + K {x — Xaiin) /2, wc cvaluatc the spring 
constant K as 



K = 



dx^ 
1 



- y^K, 



(15) 



qk = 



1 

hJi' 



with h = AttX'^H/^o. 

The condition of the vortex line formation is roughly 
^^min — 0. In such a case the pancakes will accumulate 
near the surface and form a vortex line. The magnetic 
field hi at which this happens is determined by equations 

ui = In/ii — Inui, 
In(KUi) 
1 — exp (~ui) 

which can be solved by iterations 

Ui w In hi — In In /ii , 
hi w In K + In [In (In k)] . 

Therefore, this field is only slightly larger than H^i. The 
Bean-Livingston barrier in layered superconductors with 
high anisotropy ratio is eliminated at sufficiently high 
temperatures. [22] In the field interval Hg < < iJ^i we 
may expect the accumulation of the fiuctuating PVs near 
the surface and their density would increase dramatically 
when Hz approaches H^^. 

2. Energy of pancake vortex in small-size mesas 

The previous analysis considered the case of the sam- 
ples much larger than the London penetration depth A. 
When their dimensions are comparable with A (which 
may be the case for the small mesas) the minima of the 
PVs energy near the surfaces merge and the dependence 
( 14 1 change. We may illustrate such crossover in the case 



of the strip of the width w. Using the results of the paper 
|23j where the vortex energy in the narrow superconduct- 
ing strip was calculated, we write the expression for the 



energy 
edge as 

Uoix) 



of the pancake at the distance x from the strip 



s^oiH^ii^H^^ 

47r (47rA)^ 

47T cosh(^) 



In 



W /"KX 

— sm ( — 
ttA \ w 



(16) 



In contrast to wide strips which have two energy minima 
near the edges, the energy of narrow strip at high fields 
has only one minimum at the strip center. Two minima 
merge in the center when the magnetic field reaches the 
typical value Hgo, 



HsQ 



4ti;2 



cosh 



(17) 



and above this field the minimum energy is given by its 
value at X = w/2, 



S$0 

An 



Hci - Hz 



1 - 



1 



cosh(w/2A) 



(47rA)' 



In 



At the field 



H 



<I>n 



si 



47rA2 



In 



0.5 



cosh(w/2A) 
cosh(w/2A) - 1 



(18) 



(19) 



U-min reaches zero corresponding to condition of fast stack 
formation in the center. This field is slightly larger 
than the lower critical field for the strip, H^l, which 
is determined by the energy of the complete PV stack 
at the center. In particular, for narrow strips {w < A) 
H,i « H'^l « [2$o/(™2)]{ln[zz;/«)] +0.5}, see, e.g., 
Ref. [m The formation of the PV stacks starts at the 
center when the field HgQ 
ing the approximate condition 



becomes smaller than H^l giv- 



1 



2 \ 2 



n'X 



< luK + 0.5. 



For K w 100 this happens already for rather wide strips, 
UI < 9A. Note that at this value the finite-size correction 
to the penetration field is small, ff^f k, H^^ . 



Comparing the field dependences ( |14[ ) and ( 18 ), we see 
that in the narrow strips the slope of the J7min vs 
dependence should be smaller by the geometrical factor 
[1 — 1/ cosh (w/2A)] < 1. Due to the temperature de- 
pendence of A, we may expect to observe the crossover 
between these regimes for the strips with w > A(0) by 
varying the temperature. Note also that when the min- 
imum Uo{x) just appears at the center of the strip at 
Hz = Hso, the spring constant K vanishes. 

For mesas small in both x and y directions the geo- 
metrical factor will be different, but a qualitative trend is 
expected to be the same. To illustrate this case, we con- 
sider the cylindrical geometry which also can be treated 
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analytically. The energy of the pancake vortex at dis- 
tance p from the center of mesa with radius R may be 
calculated in the same way as the energy of the vortex in 
small superconducting disk, using an image method |25j . 
In the result we have 



S<i>oH, Ig jp/X) 

47r /o (i?/A) - 



In 



R^-p" 
RX 



(20) 



where Io{z) is the modified Bessel function. The mini- 
mum energy for the fluctuating pancakes is shifted to the 
center of the mesa if the field is sufficiently large. With 
the help of expression ( 20 ) , wc find that it occurs at fields 



h,>h: 



cyl 
sO 



_$0_ 



lo 



(21) 



The lower critical field of the cylinder ' calculated 
in Ref . M 



eyl _ Iq (R/X) $0 

lo (R/X) - 1 47rA2 



luK + 0.5 - 



Kq (R/X) 
lo (R/X) 



In the limit R ^ X the lower critical field is H^f^ = 
[$o/(7i"i?^)][ln(i?/^) + 0.38] and the regime of the accu- 
mulation of the fluctuating pancakes at the center of mesa 
is always realized when the field Hz starts to approach to 
-ff^f'. For K w 100 the condition iJ^g' < H^f^ corresponds 
to the mesas with R < 5.3A. Note that numerical value 
of the critical diameter 2R « 10. 6A is rather close to the 
critical width w ~ 9X for the strip geometry. For smaller 
mesas at fields just below -ff^i' the minimum energy is 
realized for the PV located in the center, 



Uo{0) = 

■s^oHz 
47r 



47r 



{4nXy 



In 



lo (R/X) 



(22) 



Ac/A determines the Josephson length Aj = 7s. [T] The 
interaction of dense JV lattice with PV-stack lattice was 
first considered in Ref. H?) The phase distribution for 
the triangular Josephson vortex lattice at this high-field 
regime in presence of pancake vortices is determined by 
equations 



AjA0„+ sin {4>n+l-4>n + 4'v,n+l~4>v,n-kHiy-v)) 

+ sin((/)„-(/)„_i-|-0i,^„-0i,,„_i-fc/i-(y-u))=O (23) 



where kj^ — 27rsi?||/$o ^^nd (/)„^„ are the pancake- vortex 
phases. The parameter v describes displacement of the 
JV lattice. In these equations the tt phase shift between 
the phase differences in the neighboring junctions is al- 
ready taken into account. 

We consider a single pancake in the layer n = in a 
large-size sample at distance xq from the surface 



Sn arctan — 
' V a; 



y 



arctan ■ 



y 



Xq 



xo 



where the second term is the contribution from the mir- 
ror image, which is required to satisfy the condition of 
vanishing normal current at the surface d(j)^,n/dx = for 
a; = 0. At high in-plane fields the Josephson current can 
be treated as a small perturbation and Eq. ( 23 1 can be 



solved by iterations. The first iteration for the 0-th layer 
obeys the equation 



where 



A^A(/.^'^ =2cos</.,,osin[fcH(y-w)], (24) 



cos [0^,0 (r)] 



(r^ -I- Xq) — Ax'^Xq 



Without PV the first-order phase iteration is given by 
4^(r) = -2(A^A:^)-2s 
( 24 1 can be written as 



<^(r) = -2(AjA:jj)-2 sin [kniy-v)]. The solution of Eq. 



meaning that the slope of the U„ 



Hz depen- 



dence in this regime should be smaller by the factor 
[1 — l//o(i?/A)], comparing to the large mesas. This 
slope again may substantially vary with temperature due 
to the temperature dependence of A. Overall behavior is 
similar to the stripe case. For the cylindrical mesa with 
the same diameter as the stripe width the slope reduction 
is close and somewhat smaller than one for the stripe. 



C. Interaction of a pancake vortex near surface 
with dense Josephson vortex lattice 

We consider the case of the dense JV lattice when the 
JVs cores overlap. This situation realizes when the in- 
plane magnetic field -By is larger than the characteristic 
field Bcr = $0/ i'^T^Js'^), where the anisotropy ratio 7 = 



4^)(r) 



sin [kn {y - v)] 



+ 1 / rfr'ln"" 
ttXj 



ro 



(cos [(py^o (r')]-l)sin [kH{y'- v)] 

(25) 



where the integration is in the limits —00 < x',y' < 00 
and the region x' < accounts for the image contribu- 
tion. The force acting on pancake is determined by the 
phase gradient at its position. 



2s$: 



(47rA) 



c \7 
"2 V^(/)q 



(26) 



where Cq,^^ is the Levi-Civita symbol (e 



and Bo 



xyz 



-■yxz 



0) . Computing the phase gradient and taking 
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FIG. 2. (color online) Numerically co mpute d int egrals Xx(xo) 
and Ty{xo) = T'^{xq) defined by Eqs. (27c I and (27dl, which 
determine the force acting on pancake vortex near surface by 
the dense Josephson vortex lattice. 



it at the pancake position r = (xq, 0), we obtain 



^x4>\) (a;o,0) = _,2, Ix[kHXo), 



TrXjkn 
cos (kuv) 

Xq — X 



lyikHXo), 



dr " , , " ^ 



Zj,(a;o) 



(Xo - xf + y2 y^(^2 + ^2)2 _ 4^2^2 



(27a) 
(27b) 

: COS y, 

(27c) 
siny. 

(27d) 



These integrals have the foUowings asymptotics: 
2y{xo) — ?> 27r, Ix{xo) for xq — )■ and Ij,(a;o) — ?> 0, 
^xixo) — >■ 27r for — > oo. Moreover, as the force com- 
ponents have to be derivatives of the same potential, 
one can demonstrate that the integrals are connected as 
Ty{x) =T'^{x). Numerically computed integrals are plot- 
ted in Fig. ^ 

Finally, the main result of this section is the expres- 
sion for the potential energy of the PV in presence of JV 
lattice (compare with Eq. (fl])): 



U{v) = 



25$2 X^jkHx) 



cos {kHv) ■ 



(28) 



The interaction between the JV lattice and PV near the 
minimum of the surface potential is determined by the 
reduced in-plane field 

/ill = fcffA = 27rsAB||/$o- (29) 

Taking typical values A ~ 0.3/im at 70 K and i?|| = IT 
we estimate h\\ ss 1.42. 



D. Contribution of fluctuating pancakes vortices 
near surface to the conductivity of the Josephson 
vortex lattice 

Interaction with the surface pancake vortices enhances 
dissipation of the Josephson vortex lattice. Additional 
channel of dissipation gives additional contribution to the 
JV conductivity, ajv = cf'j'v + '^p^ where afy is the JV 
conductivity in the absence of PVs and dp is the excess 
conductivity due to the interaction with pancakes. The 
corresponding pancake surface contribution , gp , can be 
evaluated using energy-dissipation formula (11) 



2 W 

W — apE^LxLyS ~2—^n^Ly exp 



F2 



T 



27tT 
~1<~ 



{/2 



(30) 



H 



Here W is the energy dissipated per single superconduct- 
ing layer, and Ly are the dimensions of crystal along 
the corresponding axis and the coefficient 2 comes from 
the contribution from two yz surfaces. 

Taking in mind the relation between electric field and 
velocity for the moving JV lattice = vH^/c — 
{^o/27Tcs)uj and extracting parameters F and from Eq. 
(28), we obtain for the excess conductivity 




min 

7^ 



where 



G[f,hn 



G[T/ep,h\\ 



[l^(/l||W,nin)]^ X2(/l||M^in) 



(31) 



""min ~ In h, and T = T/sp is the dimensionless tem- 
perature. The pancake viscosity, rj, also determines the 
the flux- flow resistivity for the in-plane current, pff — 
s^()B/{c^ri). A natural scale of conductivity is the pan- 
cake flux- flow conductivity at field B — <l>o/(47rA2), 



47rA2c2ry 



Using this scale, we can rewrite Eq. (31) in somewhat 
more transparent form 



A 



4T2'^[^/£P'^II 



a 



(32) 



In spite of the rather complicated general 
ap{T,Hx:,Hz) dependence (31), the main variation 
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of CTp just below H^i is related with the change of the 
concentration of the fluctuating pancakes. In samples 
with the lateral sizes larger than A this leads to the very 
characteristic exponential CTp dependence on 



ap « A{H^j;,T) exp 



AttT 



(33) 



and the dominating temperature dependence of the pref- 
actor is given by 



These are the easiest qualitative predictions to compare 
with experiment. We remind also that in small-size mesas 
the expression in the exponent of Eq. ( 33 ) acquires a 



geometrical factor depending on mesa size, as described 
in Sec. HlBll 

A more complicated issue is the dependence of the ex- 
cess conductivity on the in-plane magnetic field. The 
general formula ( 32 ) significantly simplifies at high in- 



plane fields ^ 1. In this regime we obtain G[T, 
47rV(T^/in) and 



2^ exp 



T 



T 



3/2 \^ 



.(34) 

i.e., we expect that at high fields the excess conductivity 
rapidly decreases with increasing in-plane field oc 

We have considered the influence of the fluctuating 
pancakes on the dense moving JV lattice only near the 
yz surfaces, perpendicular to the in-plane magnetic field. 
Naturally, the fluctuating PVs also exist near the xz sur- 
faces and they also will contribute to dissipation. The 
JV lattice will displace them mainly along the x axis. 
Accurate calculations for this case are more complicated 
because one has to take into account deformation of the 
JV lattice near the surface. However, even without per- 
forming the calculation of this contribution, it is evident 
that it will be also proportional to the PVs concentration 



and therefore the dependence ( 33 ) will describe it as well. 
This dependence also should describe behavior of the c- 
axis excess conductivity for smaller in-plane fields, in the 
regime of dilute JV lattice. In this regime the fluctuating 
PVs enhance dissipation of individual JVs meaning that 
CTp cx l/Hx. 



III. EXPERIMENT 

A. Experimental setup 

Bi2212 single crystals grown by the traveling solvent 
floating zone method were fabricated into mesas sand- 
wiched by two superconducting electrodes of Bi2212 
(see Fig. ITk) with a focused ion beam (FIB) machine 



TABLE I. List of samples. L, W, and t denote length (_L H), 
width (II H), and thickness (|| c) of the mesa, respectively, see 
Figure[T^. Hp is the JV flow oscillation period experimentally 
observed. Note that Hp is slightly difi'erent from ^o/sL due 
to the ambiguity of the mesa dimensions. Units for dimen- 
sions, magnetic fields, and temperature are ^m, kOe, and K, 
respectively. 





L 


W 


t 


Hp 




Br0552 


5.0 


4.9 


0.16 


2.6 


84.7 


Br0573 


5.0 


10.2 


1.0 


2.58 


88.4 


Br0640 


1.1 


1.13 


0.25 


11.5 


84 


Br0742 


1.53 


1.45 


0.26 


8.6 


86 


Br0746 


0.8 


2.06 


0.5 


17.25 


86.5 


Br0747 


0.4 


1.5 


0.3 


34.5 


83.9 



SMI2050, SII Nanotechnology Inc. The fabrication de- 
tails are described in Ref. |28l The studied samples listed 
in Table [J include two medium-size mesas with lateral 
sizes > 5 fim and four small mesas with lateral sizes ^ 1 
/im. 

Resistance of the mesa was measured with the four- 
probe method with either dc or ac current excitation. 
For the ac excitation, the lock-in detection was used for 
voltage measurements. External magnetic fields H were 
applied by a split-pair superconducting magnet that gen- 
erate a horizontal magnetic field up to 80 kOe and the 
angle between the magnetic field and the a6-plane (a;- 
axis) was varied by rotating the sample probe with a 
precision rotator. Perpendicular (c-axis) component of 
the magnetic field is given by HsinO. In all plots 
presented here, is always swept from negative to pos- 
itive, meaning that PVs are expelled from and introduced 
into the mesas in the negative and positive regions, 
respectively. Therefore asymmetry with respect to sign 
of Hz in some plots may be caused by hysteretic behav- 
ior. 



B. Experimental results and discussion 

The c-axis resistivity pc in the vicinity of 9 = gives 
the JV flow resistance. With increasing H^, Pc suddenly 
drops due to penetrating PVs. Without taking fluctuat- 
ing PVs into consideration, the value of the field when 
Pc drops should correspond to the lower critical field Hd 
(corrected by the demagnetization factor). The demag- 
netization factor of a small mesas (L or VF ~ l^m) is 
small because the thickness of the whole crystal is about 
10 /im for all samples and it is larger than the width. 
For the midsize mesas (L or VF ~ 5 — 10/im), due to 
demagnetization effect, the field at the edge of junctions 
in the Meissner state should exceed the external field by 
factor ^1.2-1.3. An example of pc as a function of at 
different applied fields for one of the medium-size mesas 
is plotted in Fig. [3](a). We see that the steep drop in pc 
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FIG. 3. Typical dependences of the resistivity pc{Hz) in the 
midsize mesa Br0573 at 60 K (a) and small-size mesa Br0746 
at 70 K (b). (a) Sharp drop corresponds to penetration of 
the PV stacks. We focus on behavior below this drop, (b) In 
small-size mesas, in addition to the initial sharp drop, quasi- 
periodic drops due to penetration of the individual PV stacks 
are found. In the case of either small magnitude of tilting 
magnetic field (2 kOe) or magnetic field parallel to the c-axis 
(solid line), stepwise changes in pc indicated by arrows are 
observed. 

in the low-field region is smeared with increasing parallel 
field H. We focus here on the quantitative characteriza- 
tion of behavior slightly below this drop. 

Field behavior of small mesas has two important fea- 
tures. Firstly, when the magnetic field is aligned with 
the layers, due to strong interaction with the edges, 
the dense JV lattice has rectangular structure almost at 
all fields and only transforms into the triangular con- 
figuration in the vicinity of fields corresponding to in- 
teger number of fluxes per junctions, Hn — nHp with 
Hp ~ <i>o/sL.[THl[in] Values of Hp for different mesas are 
listed in the Table|l] The c-axis resistivity has pronounced 
oscillations as function of Hj. with maxima located near 
ffn-Hmil] Secondly, when the c-axis magnetic field is ap- 
plied, penetration of the individual PV stacks can be re- 
solved. In the case of magnetic field tilting from from the 
a6-plane, a PV stacks penetrating into the mesa strongly 
impede JV flow along the afo-plane, which causes a sharp 



FIG. 4. The dependences of the excess conductivity on Hz for 
the sample Br0573 at T = 60 K in (a) 2 < < 10 kOe and 
(b) 15 < H < 60 kOe. Solid lines describe exponential depen- 
dence ( |33[ ). One can see that the data follow the theoretical 
prediction. 



drops in pc. With increasing Hz, the separation between 
the drops decreases and finally corresponds to the period 
of one flux quantum per mesa area, ^q/LW , as shown in 
Figj3](b). These quasi-periodic drops in pc are naturally 
attributed to penetrations of the individual PV stacks 
into the mesa. Even in the absence of JVs {H \\ c), 
stepwise increase in p^ (shown as a solid line) pointed 
by arrows in Fig. [s] (b) with similar periodicity was 
observed. These steps are attributed to additional in- 
terlayer phase fluctuation caused by the penetrated PV 
stacks. [28 Therefore in small mesas the c-axis resistivity 
has oscillating dependences on both field components. 

In order to quantify the additional dissipation of 
the JV lattice due to penetration of fiuctuational PVs, 
we show the plots of the excess c-axis conductivity 
ap{H^, Hz) = ac{Hx, Hz)- <Jc{Hx,0) a.s a function of the 
c-axis magnetic field Hz for several samples at different 
temperatures in the following figures. Here, ac{Hx,0) = 
1/ Pc{Hx,0) is the conductivity for the magnetic field 
aligned with the layers, at = {Hz ~ 0). Only the data 
at temperatures higher than T — 60K are shown because 
at low temperatures (below 50 K), the surface barrier 
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and bulk pinning for PVs cause pronounced hysteresis in 
the dependences pc{Hz), which makes an analysis quite 
complicated. 

Figure |4] illustrates the evolution of the dependence 
of the excess conductivity <7p with increasing field ob- 
tained from data presented in Fig. [3] Other mesas show 
similar behavior in the penetration region. In the figure, 
the predicted exponential dependence ( |33| is plotted as a 
thick solid line for each T. Since the constant A is the sole 
fitting parameter, the slope of the semi-log plot of the ex- 
perimental data provides a direct test for Eq. ( 33 1 . The 



first general observation is that the experimental data 
clearly show the existence of the exponential regime just 
before the penetration of the perpendicular field into the 
mesas. This exponential increase spans up to five orders 
of magnitude. The slopes of the ln((Tp) vs curves are 
in good agreement with Eq. ( 33 1 . The observation of 



this regime provides strong evidence of generation of the 
fluctuating PV near the sample surface. 

There are several common features in the experimen- 
tal data. The slopes of dependences In dp vs remain 
independent of the in-plane field only up to the typi- 
cal fields ~40 kOe, in agreement with theory. At higher 
fields (not shown) the slope starts to decrease and at the 
same time the experimental dependencies ap{Hz) are no 
more symmetric over ±Hz ■ Probably at such a high mag- 
netic field it is difficult to reliably control the small 
component. At small fields < 10 kOe the excess conduc- 
tivity is only weakly depends on the field while at higher 
fields it starts to decrease rapidly. This behavior is most 
probably related with the crossover into the dense-lattice 
regime. Figure [5] shows dependence of the preexponen- 
tial factor A on the in-plane field for mesas Br0552 at 
T = 70K and Br0573 at 60 K. One can see that decay at 
fields Hx > 6kOe is close to theoretical dependence 
given by Eq. (34). 



1. Small mesas: triangular vs rectangular JV lattice 

An important feature of small-size mesas is the exis- 
tence of regions of the rectangular vortex lattice at fields 
Hx > BcrL/ (7s). In this range with increasing magnetic 
field the JV lattice undergoes a series of structural phase 
transitions between the triangular and rectangular con- 
figurations, with triangular structures located around the 
fields nHp and rectangular structure located around the 
fields {n + l/2)Hp. The widths of the triangular regions 
shrink with increasing field |181 119j . The triangular con- 
figuration is expected to have much stronger interaction 
with the PVs than the rectangular lattice because there 
are no in-plane currents in the latter state. 

Figure [6] presents the dependences of the excess 
conductivity for two small-size mesas Br0746 and Br0742. 
The panel (a) shows data for mesa Br0746 at 75 K for 
fields smaller than Hp — 17.26 kOe. We can see that 
there is again a wide range of exponential dependence and 
the slope dlnap/dHz agrees with theoretical predictions. 
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FIG. 5. The normalized preexponential factor A in Eq. 1 33 1 
as a function of in-plane field H^c for the mesa Br0552 at 70K 
and Br0573 at 60 K. The factor is normalized by its value at 
H^=20 kOe. 



Figure [6]jb) shows the excess conductivity for the small 
mesa Br0742 for three values of field 4.3, 8.6, and 12.8 
kOe almost corresponding to 0.5, 1.0, and l.hHp. At 
the last field the rectangular lattice is realized. We can 
see that while agrement with the predicted dependence 
is very good for the first two values of field, for the last 
field the excess conductivity before penetration is much 
smaller which makes the transition much sharper. This 
confirms that the fluctuating pancakes have much weaker 
influence on dissipation of the rectangular vortex lattice. 

We finish this subsection with a small comment on the 
oscillating behavior of the c-axis resistivity in small-size 
mesas at larger due to the penetrating individual PV 
stacks illustrated in Fig.jsJ^b). This behavior actually ap- 
pears due to the competition between two opposite fac- 
tors. On the one hand, the PVs in the stacks directly con- 
tribute to the dissipation of the JV lattice reducing the 
resistivity. On the other hand, penetration of the mag- 
netic flux inside the mesa relaxes the Meissner currents 
near its surface which reduces the contribution coming 
from the fluctuating PVs and increases the resistivity. 
We can see that slightly above the penetration field the 
resistivity increases meaning that the second mechanism 
actually prevails in this region. Moreover, continuous in- 
crease of the resistivity for H = IQ and 14 kOe indicates 
that the full PV stack in the mesa center is not formed 
at once but continuously builds within finite field range. 



2. Evolution with increasing temperature and finite-size 
effect m small mesas 

Figure [7] presents evolution of the dp vs Hz depen- 
dences with increasing temperature at fixed field for two 
small-size mesas. In both cases fields are close to the 
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FIG. 6. The excess pancake conductivity for two small-size 
mesas, (a) Data for the mesa Br0746 at 75 K for three values 
of magnetic filed smaller than Hp = 17.25kOe. Theoreti- 
cal exponential dependences and corresponding preexponen- 
tial factors are also shown, (b) The excess conductivity for 
the sample Br0742 at 70 K also for three values of field. For 
the first two values 4.3 and 8.6 kOe corresponding to Q.^Hp 
and l.QHp agreement with theory is very good. There is no 
agreement for the last value, 12.8 kOe ~ 1.5-ffp, the transi- 
tion is much sharper and occurs at higher H^- At this value 
the rectangular JV lattice is realized in the mesa which is less 
sensitive to the fluctuating PVs. 



FIG. 7. Temperature evolution of the excess conductivity 
for small-size mesas, (a) Excess-conductivity plots for the 
small-size mesa Br0640 for the field 11.8 kOe corresponding 
to the triangular lattice at four values of temperature, 60, 70, 
75, and 80 K. Solid lines represent exponential dependences 



given by Eq. (33 1 at temperatures for corresponding colors. 



whereas the broken lines are given by Eq. ( 35 I for the small- 
size regime. One can see that the exponential dependencies 
agree with the prediction for the large-size regime at 60, 70, 
and 75 K but the slope is considerably smaller at 80K. (b) 
Data for the mesa Br0747 from 60 to 82 K at 3 T. The ex- 



ponential dependence for the small-size regime ( 35 1 shown by 
broken lines describes data very well. 



corresponding values of Hp. We can see that for the 
mesa Br0640 (panel (a)) these dependences agree with 
the prediction for the large-size regime for temperatures 
60, 70, and 75 K but slope dlnap/dHz becomes consider- 
ably smaller for 80K. For the smallest mesa Br0747 with 
the width L = 0.4 /im (panel (b)) the slopes are smaller 
than predicted by Eq. ( 33 ) for all temperatures and de- 
crease with increasing temperature. This behavior can 
be naturally explained by the finite-size effect considered 
in Sec. Ill B 21 As the other size of this mesa W is 3.75 
times larger than the width L, we can apply the strip ap- 
proximation. In this case when the London penetration 
depth A becomes comparable with the mesa smallest lat- 
eral size, the energy minimum of the PV, J/mim is shifted 
to its center and the slope the JJminVS decreases, thus 



Hz must be substituted by H^ [l — cosh ^ (L/2A)] : 
(Tpf^ A{Hx,T) exp I — . 1 — cosh" 



47rr 



(35) 



This means that from the slopes we can extract the 
temperature-dependent A. Figure [S] (a) presents the 
temperature dependence of obtained in such a 

way. It shows expected linear dependence near Tc, 
A-2 = - T/T^P), with the Ginzburg-Landau 

value Aq — 0.1/im and the mean- field transition tem- 
perature « 86K which is 2 K higher than the 
resistivity transition temperature for this mesa, Tc = 
84K. Additional information can be obtained from the 
temperature dependence of the preexponential factor 
A presented in the right part of Fig. [8] As follows 
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FIG. 8. (a) Temperature dependence of extracted from 
slopes dXnojUdHz for mesa Br0747 using the activation en- 
ergy, Eq. (181. (b) Arrhenius plot of the preexponential fac- 
tor A for the same mesa. The origin of huge value of Aq is 
explained in the text. 



However for this mesa the sizes W and L are close so 
that the narrow strip model is not applicable. Accurate 
description of this case requires numerical analysis. 
IV. CONCLUSIONS 

Layered superconductors in magnetic fields are charac- 
terized by complex interplay between the Josephson and 
pancake vortices. In particular, dissipation of the Joseph- 
son vortex lattice is proven to be extremely sensitive to 
the presence of pancake vortices. We utilize this property 
to probe the fluctuating pancake vortices. In the Meiss- 
ner state, in the c-axis magnetic field smaller than the 
lower critical field, the pancakes vortices can not form 
the Abrikosov vortex lines. In the fiuctuational regime, 
however, the individual pancake vortices may exist near 
the surface or inside the small samples. They lead to 
additional contribution to the c-axis conductivity of the 
Josephson vortex lattice which has very characteristic ex- 
ponential dependence on the c-axis magnetic field. While 
in mesas with lateral sizes significantly larger than the 
London penetration depth A the slope dhi{ap)/dHz of 
this dependence is universal, in smaller mesas it acquires 
the geometrical factor depending on the ratio size/A. We 
systematically studied the excess c-axis conductivity of 
the JV lattice due to the fiuctuating pancake vortices in 
mesas fabricated out of Bi2Sr2CaCu208+5 crystals. The 
predicted exponential dependence is clearly observed and 
its slope agrees very well with the theoretical value. Ana- 
lyzing the temperature evolution of the slope in the mesa 
with smallest width and extracting the geometrical fac- 
tor, we were able to restore the temperature dependence 
of the London penetration depth. These findings pro- 
vide strong evidence for the existence of the fiuctuating 
pancakes in the Meissner state and demonstrate that the 
Josephson-vortex lattice provides a unique tool to probe 
them. 



from Eq. (fTsl), it has Arrhenius temperature depen- 
dence, A =~3o exp(— [/o/^bT")) with the activation en- 
ergy Uo = [s<I>g/(47rAo)2][ln(L/7rf) + 0.5]. This agrees 
with experiment and the fit gives Uo/ks = 1.08 • 10'*K. 
This is consistent with the above value of Aq if we as- 
sume ^ w 50nm at 75K. Note also that the huge value 
of Aq obtained from the fit, Aq « exp(125)[ricm]~^, is, 
in fact, very reasonable because Aq contains the fac- 
tor exp{Uo/kBT^^) « exp(126). The obtained value of 
Aq for our overdoped small mesa, however, is somewhat 
smaller than the values reported in the literature. [21] 
Note that the mesa Br0640 also shows the same trend. 
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